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Abstract 

We analyze the properties of mesons in 1+1 dimensional QCD with bosonic and 
fermionic "quarks" in the large N c limit. We study the spectrum in detail and show that 
it is impossible to obtain massless mesons including boson constituents in this model. We 
quantitatively show how the QCD mass inequality is realized in two dimensional QCD. 
We find that the mass inequality is close to being an equality even when the quarks are 
light. Methods for obtaining the properties of "mesons" formed from boson and/or fermion 
constituents are formulated in an explicit manner convenient for further study. We also 
analyze how the physical properties of the mesons such as confinement and asymptotic 
freedom are realized. 



email: a00500@cc.hc.keio.ac.jp 
2 email: tomo@th . phys . titech .ac.jp 



1. Introduction 



In confining gauge theories, physically observable particles at low energies have no 
gauge charge and are bound states of the charged matter that appear in the gauge theory 
Lagrangian. A prototype of such a theory is QCD wherein the gauge group is SU(3) and 
the charged matter fields are vector fermions in the fundamental representation. Several 
natural generalizations of QCD exist; we may use gauge groups other than SU(3), we may 
use chiral fermions or we may choose representations more complicated than the funda- 
mental representation for the matter fields. In this case, we cannot put in too many or too 
large representations if we want to preserve asymptotic freedom. Another generalization 
is to consider boson constituents as well as fermion constituents. 

It is this last, as well as the first generalization in two dimensions in the large iV c limit 
that we shall investigate in this work. The possibility of boson constituents arise necessarily 
in several contexts, such as technicolor with multiple stages of symmetry breaking Jl[], 
QCD or technicolor with super symmetry 0, so called bosonic technicolor as well as 
the Standard ModelQ. When both boson and fermion constituents exist in the theory, 
"meson" states of both Bose and Fermi statistics may arise. In general, it is difficult to 
derive the properties of the bound states from first principles. By using the large iV c limit 
in two dimensions, we may analyze the properties of these meson states concretely. 

Two dimensional QCD in the large N c limit has greatly contributed to our current 
understanding of the gauge theory dynamics by providing us with a model where the 
properties are explicitly calculable analytically. Also, two dimensional QCD has been used 
to test the validity of various approaches and approximation schemes applied to QCD. 
The model was first solved by 't Hooft |5| and some further physical properties such as 
some current matrix elements, the asymptotic freedom of mesons were studied in some 
subsequent works ||0. The formulation was extended to include boson-boson bound 



states HH and boson-fermion bound sates 110]. Mesons made only from fermionic quarks 



or the bosonic quarks obey Bose statistics but the boson-fermion bound state obeys Fermi 
statistics. (Hereafter, often referred to as ff, bb and bf cases.) 

In this work, we shall extend the investigation of the physical properties of the mesons 
made from fermions and generalize the results to the mesons made from bosons only and 
bosons and fermions. The results will enable us to compare the three cases and see the 
differences and similarities that arise between mesons made from constituents of various 
statistics. The spectrum of mesons is investigated both analytically and numerically. We 
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will study the case when the quarks are heavy analytically. For mesons involving light 
quarks, we establish a number of results analytically and further analyze the problem 
numerically. The numerical approach will be formulated explicitly in all the three ff, bf 
and bb cases which should be useful for further study. 



We use these results to see how the QCD mass inequality]! 1 ] |12|] applies to two 
dimensional QCD in the large iV c limit. Given two different types of quarks, the QCD 
mass inequality states that the meson made from the same quarks is on average lighter 
than the meson made from different quarks. This non-trivial inequality, however, does not 
tell us by how much these meson masses differ, a question which we are able to answer 
analytically in some cases and numerically for all quark masses. Also, while there is 
no reason to doubt this important inequality, as was pointed out in the original articles 
themselves, the inequality is not completely rigorous. We find it satisfying to be able to 
study how the inequality is realized in a simplified version of QCD. To our knowledge, 
the QCD mass inequality has not been previously shown to be satisfied in such a concrete 
manner. Here, the inequality is applied to ff, bf and bb cases. Except in the bf case, 
the inequality is necessarily an equality when the constituent quark masses are the same. 
In the bf case, such needs not be the case and indeed we find that it is always a strict 
inequality in the bf case. 

The paper is organized as follows. In section 2, we briefly review the integral equation 
satisfied by the wave function of the meson in the large N c limit, partly to fix the notation. 
In section 3, we analytically study the static properties of mesons when the quark masses 
are large. We further formulate the bound state equation as a linear eigenvalue problem 
for general values of the quark masses and analyze the problem numerically. In sect 4, we 
analyze some matrix elements and see how confinement and asymptotic freedom is realized 
in the bf case. We conclude with a brief discussion. 



2. Wave function of mesons 

In this section, we briefly summarize the equations meson wave functions satisfy and 
some basic properties of the solutions in the ff[[|, bf|L(J and bb0 cases. The classical 
Lagrangian of QCD coupled to fermions and bosons is 



Both fermions and bosons are in the fundamental representation. We shall refer to the 
fields i/jf,(f>b as "quarks". We fix the gauge to the light-cone gauge A_ = A + = 0, where 
(a 1 ±a°)/V2. Light -cone gauge has the advantage that there are no gluon self 
interactions in 1+1 dimensions. 

We take the large N c limit by letting N c go to infinity keeping g 2 N c and the quark 
masses to be of 0(1). There are contributions of 0(1) to the quark propagators which 
compete with the classical contribution. These contributions may be incorporated by 
solving the Schwinger-Dyson equations recursively to obtain the full propagators as follows: 
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We note that the quantum corrections to the mass is identical both for the fermionic 
and the bosonic quarks. In this work, we will use an infrared cutoff A_ which seems to 
be more convenient for deriving the physical properties of the mesons. Other infrared 
regularizations may be used, but of course do not affect physical results || 

A meson is formed as a quark anti-quark bound state and its wave function satisfies 
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Here, = m 2 h ix/(g 2 N c ), namely the mass squared of the quarks counted in the units 
of the QCD scale. These equations may be obtained by summing the graphs of the ladder 
type. Here and below, we refer to the solution of the Bethe-Salpeter equation, v?(x), as the 
wave function of the meson, in analogy with the non-relativistic case and in accordance 
with the previous literature, x = p_/r_ is the momentum fraction carried by the quark 
and 1 — x by the anti-quark in the infinite momentum frame. In the 1/N C expansion, sea 
quarks are suppressed by 1/N C and gluons have no physical degrees of freedom so that the 



all the momentum of the meson is carried by the quark and the anti-quark. Pf denotes 
the principal part integral defined by 



pjdx f(x) = — Jdx 



f(x + ie) + f(x - ie) 



(2.5) 



In the bf and the bb cases, the above meson equations are not Hermitean with respect to 
the standard measure Ldx and we shall often use the conjugated equation 
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The appropriate conjugation is 
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A formula we find useful is 
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Note that the expression no longer involves the principal part integral. From this equation, 
it immediately follows that in the ff case when the bare quark masses are positive (f3 a , fib > 
0), the meson mass is positive [[5]. In the bf, bb cases, it is clear that the meson mass is 
positive when flat fib > 1- Another property that may be obtained from this formula is that 
for the same quark-antiquark masses, the mass of the meson is the lightest for the bb case, 
heaviest in the ff case. Here and below, we often refer to the bare masses of the quarks 
m/ ; b as quark masses. It should perhaps be commented here that the quark masses m/^, 
or equivalently are not directly physically observable and that they receive quantum 



corrections. However, we note that the quantum corrections to both the fermionic and the 
bosonic quark are identical so that it is consistent to compare the quark masses on the 
same footing. 

The wave function vanishes at the boundary as (f(x) ~ x 7/i where x denotes the 
momentum fraction of the fermionic or the bosonic constituent. 7/^ is determined by the 
equations 

f (3 f - 1 + tt 7/ cot(7r 7/ ) = f , 2n . 
\ (3 b - 1 - itj b tan(7T7 6 ) = b 

Since the meson equations are Hermitean, the spectrum is real. The widths of the 
mesons are of 0(iVfl avors /iV c ) and are suppressed in the large N c limit. For higher mass 
states, the approximate wave functions and the respective masses are 

(p(x) ~ sin7r/cx, (meson mass) 2 ~ g 2 N c nk, k = 1,2,3... (2-12) 

resulting in a linear Regge type trajectory for large k. Here, k labels the fc-th lightest 
meson bound state formed from the quark-antiquark pair. We note that the relativistic 
effects are important in this problem, even qualitatively; in the non-relativistic case, the 
bound state mass squared behaves as k 2 ^ 3 rather than linearly. 



3. Static properties of mesons 

In this section, we obtain and analyze the static properties of mesons using both 
analytical and numerical methods. 

3.1. Heavy quarks 

When the quark masses are heavy, we may use the variational method to obtain the 
meson wave function analytically. This method was previously applied to the ff case when 
the quarks and anti-quarks have equal mass |J. Here, we extend the analysis to include 
the case when the quark and the anti-quark have different masses and further generalize 
this to the bf and bb cases. We refer to the "heavy quark" when the quark mass is large 
compared to the QCD scale, (m ^> g <^ /3 > 1). 

Using the trial function 

<p^( X ) = ^y /4 e -^-o?/2 (31) 
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we obtain 



(V, Hip) = ^ + + ( % + ) 1 + V^c + J(x ) + 0(Pc~ 2 ) (3.2) 

Varying with respect to xq and c we obtain 

C = (TT/3^)- 1 / 3 (yp~ a + V 7 ^) 873 (1 + 0{C~ 1 )) 
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Here, we have obtained a clear physics picture of the meson as a bound state of the quark 
and anti-quark with the momentum distributed in proportion to their masses. As the q, q 
masses become larger, the momentum distribution becomes narrower as /3 -1 / 3 and this 
simple constituent quark picture becomes more accurate. 

Since we are considering different masses for the quark and the anti-quark, it may 
seem we should analyze the effect of including a variational function that is asymmetric 
with respect to the center of the wave function, such as (x — xo) exp — c(x — xq) 2 /2. It can 
be immediately shown that to the order we considered above, including this function in 
the variational problem does not change the results at all. 

To systematically obtain the wave function when the quark masses are heavy we may 
use an orthonormal basis for the wave functions 
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where H n is the n-th Hermite polynomial defined by 
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To determine xq, c to leading order, as in ( |3.3| ), we need the matrix elements of the Bethe- 
Salpeter equation to C(/3 1 / 3 ). To this order, the statistics is unimportant. In all the three 
ff, bf and bb cases, the matrix elements of H are 
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The one dimensional variational problem just using ip^ j q {x) is not quite enough to determine 
the lowest meson mass to the order in /3 -2 / 3 given above in ( [3.4|) . It turns out, however, that 
only when we consider trial functions (p^ q (x), with n > 4 that the meson mass is affected. 
When we include the variational functions with n > 4, we find that the physical effect is 
to decrease the numerical coefficient of the second term in the square brackets of ( |3.4| ) by 
less than 1%. The last term, w s , is not affected at all. The smallness of the corrections 
is not surprising; while the contribution of more complicated variational functions are not 
suppressed by powers of /3 -2 / 3 , they are suppressed numerically just as the contribution 
from the higher excited states to the ground state energy are suppressed in perturbation 
theory in quantum mechanics. This suppression is strong since there is a contribution 
to the meson mass only when we include variational functions with n > 4. A couple of 
comments are in order; first, the boundary conditions at x = 0, 1 are not exactly satisfied in 
this approach, so that the expansion is not completely systematic. However, since the value 
of the wave function at the boundaries is exponentially suppressed as exp — const, x /? 2 / 3 , 
the approximation should be good when the quarks are heavy. Second, the growth of 
the matrix elements with the increase in the dimension of the variational space as seen in 
( p.8| ) indicates that the expansion in /3 -2 / 3 is an asymptotic series. This situation is quite 
common in quantum theories [0 . Further analyzing the problem through this variational 
approach does not change the physics picture obtained above and we shall not pursue this 
further here. 



3.2. Variational methods using powers of the momentum fraction 

An effective method, both anaiytically and numerically, for obtaining approximate 
solutions to the integral equation for the meson wave function, (|2 . 3| ) , ( pT6|) is to reduce 
the problem to a finite dimensional diagonalization problem. An example of a variational 
approach particularly effective in the case of heavy quarks was given above. When the 
quark masses are light, a variational scheme that is effective is to use a basis 

v 2k = [x(l-x)] 7+ \ v 2k +i = [x(l-x)] 7+fc (l-2x), k = 0,1,2,... (3.9) 
This scheme was previously used in the ff case|I3] and more recently in the investigations 



of the Schwinger model [TIJ. Below, we treat the ff, bb cases when the quark and the 



anti-quark masses are identical ((3 a = (3b = (3). In this case, the integral equation ( |2"1)| ) is 
reduced to finding the solution of the generalized eigenvalue problem 

det (H ZJ - ii 2 N l3 ) = 0, i, j = 0, 1, 2, . . . (3.10) 

where 

Hij = (vi, Hvj) = Kij + Uij, Nij = (vi, vj) (3.11) 
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K 2k , 2l = (P-l)B(2j+k+l, 2 7 +fc+Z), K 2k+1;2l+1 = 47 + i 2 k ~+li + 1 B ^ +k+l ^ 2 ^+ k + l ) 

(3.13) 

2(2 7 fc + ( ig/) ) g (7 + 7 + k)B{l + 1,1 + 1) ff 

8kl+(8l t( 2 f + t2n 8l2+2l B ^ + k + 1/2, 7 + k + 1/2)S( 7 + I + 1/2, 7 + I + 1/2) bb 

(3.14) 

" +W I ^gg!? 7 % + k + 1/2, 7 + k + l/2)fl( 7 + / + 1/2, 7 + ^ + V2) bb 

(3.15) 

H^, — when i + j = 1 mod 2. The boundary conditions mentioned in the previous 
sections require that 7 > which guarantees that the meson masses and the matrix 
elements given above are finite. 
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It is illuminating to study the simplest one dimensional case analytically. In this 
case, the meson mass squared is 

4B(2 7 +1,2 7 +1) tt 

U+l.) ^(7+1/2,7+1/2) bb 
W + 4^ B(2 7 +l,2 7 +l) DD 

We may vary 7(> 0) to obtain an upper bound on the meson mass. In the ff case, when the 
quark mass is positive ((3 > 0) the problem is well behaved and it was already established 
that the meson mass is positive using ( |2.9| ). When (3 = the meson mass is zero. When 
(3 < 0, (H) is not only negative but is unbounded from below. The situation parallels that 
of the 1/r 2 potential in quantum mechanics; when the coupling is smaller than the critical 
value, the problem becomes ill behaved. In the bb case, this critical coupling is at j3 = 1; 
when (3 < 1, the meson mass is unbounded from below and is not acceptable physically. 
As a result, the bosonic quark mass needs to satisfy the condition (3 > 1, which we shall 
adopt from now on. When (3 > 1 the meson mass is positive and the problem is well 
behaved. Using this variational approach, we may establish that < [i 2 < n 2 /4 for (3 = 1. 
Consequently, it is not possible to obtain a massless bb meson in this model. Combined 
with the QCD mass inequality to be explained below, this excludes the possibility of 
obtaining massless mesons with boson constituents. 

We may use this formulation for the numerical computation for the spectrum. In the 
bf case, it is natural to choose 7 to satisfy the boundary behavior (|2.11|) . This formulation 
is most convenient for the ff case with light quark masses (3 < 1. 

The basis we have chosen is not orthonormal. An orthonormal basis is 

v[ = z**>(l - x)^G n (2 l0 + 2 7l + 1, 2 7o + 1; x) (3.17) 

where G n (a, b; x) denotes the n-th Jacobi polynomial. While we may compute the matrix 
elements of the diagonalization problem analytically using this orthogonal basis, we have 
not found a compact general formula for the matrix elements as in (|3.10|) and we shall not 
use this approach here. 

3.3. Multhopp's method 

Another numerical method we shall employ, which is sometimes called Multhopp's 
wing dynamics, is valid for all quark masses in the three ff, bf and bb cases. This method 
has previously been applied to the ff case with mesons formed from quarks with identical 
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masses m 



where 



16 . Here, we apply this method to a meson made from a quark and an anti- 
quark of different mass in the ff, bf and bb cases. We approximate the meson wave function 
by 

<fi(x) = a n sin n6>, x = — - (3.18) 
n=l 

Using integration by parts, the integral equation ( |2.3j ) for the meson wave function reduces 
to 

K 

fi 2 ^ = Y,M n (9)a n (3.19) 

n=l 

M nW S 2 ( + Alij) Si nn0 + 2n ( 2*2£» + B n (e)) (3.20) 
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ff 

. cos $ cos n6 

B n (0) = <( 2(i+cose) bf (3.21) 

-cos cos n#+<5„.i /8 
sin 2 9 

We solve the equation by evaluating it at the points 6 = 9 n: n = 1, 2, . . . , K where we 
used the notation 9j = ixj / (K + 1). Using the relation 

K K + l 

^ sin 9 nk sin 9 mk = — - — 5 nm (3.22) 

k=i 

the problem reduces to a finite dimensional matrix eigenvalue problem 

2 K 

^2 A nk<ik = (J> 2 a n , A mn = — — sin 9 mk M n (9 k ) (3.23) 
k k=i 

We now comment on the relative merits of the numerical schemes we used. First, we 
note that when solving the problem numerically, it is beneficial both in terms efficiency 
and in terms of accuracy to be able to compute the matrix elements of the diagonalization 



problem analytically. If one does not insist on this property, we may use (|2.9[ ) to compute 
the matrix elements by numerically performing the double integral for any basis satisfying 
the boundary conditions. Multhopp's method is robust in that it tends to produce matrix 
elements of order one even when K is large and does not give rise to almost singular 
matrices numerically. Also, the diagonalization problem is well behaved for all values of 
the quark masses. Except in the region of light quark mass (f3 < 1) the convergence is 
sufficiently fast; K = 200 is more than enough to obtain the meson mass to eight digits. 
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K = 800 is a reasonable computational task requiring the order of an hour of cpu time 
on a current workstation. For the ff case when the quarks are light and are of equal 
mass, the basis involving powers of the momentum fraction is useful. A ten dimensional 
diagonalization suffices to obtain the spectrum to five significant digits. However this 
approach tends to lead to almost singular matrices numerically as we increase the dimension 
of the variational space. It may be possible to overcome this problem by a clever choice of 
basis. The convergence is slower than the Multhopp's method except in the ff case with 
light quark masses of equal mass. 

3.4- The meson spectrum and the QCD mass inequality 

In all the three ff, bf and bb cases, the spectrum leads to a linear Regge trajectory 
for the higher mass states as we saw in Q2.12| ). Here, we plot the examples for three sets 
of parameters, (i) (3 a = fa = 1, (ii) (3 a = 1, fa = 10 and (iii) (3 a = /3& = 10 for ff, bf and bb 
cases in fig. 1. Here and below, unless otherwise stated, the numerical errors are too small 
to be visible on the plots. 




0+ — * — * — 1 — 1 — 1 — 1 — * — 1 — 1 — 1 — 1 — * — * — 1 — 1 
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 

fig. 1 The spectrum of mesons for the cases (i) (3 a = flb = 1, (ii) fta = 1, A> = 10 and 
(iii) P a = [3}, = 10. The horizontal axis is labeled by the number when we count from the 
lightest meson ( "principal quantum number") and the vertical axis is labeled by the meson 
mass squared counted in the units of the QCD scale (g N c /tt). The spectrum for the ff, bf 
and bb cases are joined by solid, dashed and dot dashed lines, respectively. The ff, bf and 
bb cases are barely distinguishable for (i) and indistinguishable in the other cases. 

Let us now analyze how the QCD inequality is realized in QCD in 1 + 1 dimensions in 
the large N c limit. The QCD mass inequality was originally shown for the four dimensional 
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case 



I 1 j [12 . At a formal level, the inequality may be extended to the ff, bb and bf cases 
in two dimensions. Given the analytical form of the masses when the quark masses are 
heavy, we may analyze how the QCD mass inequality is satisfied in this case. Comparing 
the masses, we obtain 



2/iab — (fJ-aa + Hbb) 



3^ j"/ ^ ^1/3 



(^/K+v%) - 2 2/3 (pl /6 + A 



1/6 



+ 0(/3~ 1 / 2 ) > 
(3.24) 

to this order in all three ff, bf and bb cases. We see that the mass difference arises not at 
leading order but at next order since the leading order only contains the rest masses of the 
quarks. Furthermore, we see that the mass inequality can be an equality to this order only 
if the masses of the quarks are the same. In the ff and bb cases, this inequality is an exact 
equality when the masses of the bound quark and anti-quark are the same ((3 a = (3b) since 
the mesons being compared can be the same meson. When applying the mass inequality 
to the bf case, this is not the case so that it can be and is a proper inequality. In the bf 
case fi a b in the inequality is the mass of bf fermionic meson and the \i aa , \ibb are the masses 
of the ff, bb mesons. When (3 a = (3b(= (3), we may compute to next order and obtain 

1 / r- 7T\ 0.022 

1\lab - {Vaa + Vbb) = -^(V2-l--)~ —jp- > (3.25) 

As mentioned above, since the expansion in (3~ 2 ^ is an asymptotic one, this analytic 
derivation of the inequality is not rigorous and perhaps should be considered illustrative. 
The numerical results below clearly show that the inequality in the bf case is a proper 
inequality. 

The QCD mass inequality may be analyzed numerically for arbitrary values of the 
quark mass. We plot the normalized mass inequality ([i a b — (^ aa + A*6b)/2) / fJ-ab, which is 
dimensionless, against (3b in fig. 2, fig. 3 for (3 a = 1, 10 respectively. We immediately see 
that the QCD mass inequality is a proper inequality in the bf case. The values for the 
normalized mass inequality at f3 = 0, while not visible on the plots, are finite and can go up 
to 0.3 and 0.1 for the (3 = 1 and (3 = 10 cases, respectively. Also, even though the relative 
mass inequality is increasing as we increase (3b, at some point, the relative inequality starts 
to decrease and is never more than 0.012, 0.002 for (3b > (3 a in the (3 = 1, 10 cases respec- 
tively. For (3 = 1 the numerical data can have appreciable errors as shown on the plots. 
The errors were crudely estimated as follows: Since the numerical data has not completely 
converged, we linearly extrapolated the finite dimensional results by using the analytically 
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known (3 a = (3b = ff case as a guide. The error in the extrapolation was estimated by the 
statistical error in the extrapolation when various sets of data were taken. 
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fig. 2, fig. 3 The relative QCD mass inequality for the cases (3 a =1 and 10. fb (bf) means 
that (5 a is the mass of the fermionic (bosonic) quark. 

One property that strikes the eye is that the QCD mass inequality is surprisingly 
small, even when the quark masses are of order one in QCD scale. Here, we remind the 
reader that the QCD scale is taken to be g 2 N c /ir, which is a conservative choice when 
we consider that the Regge slope parameter is larger by a factor of n 2 as in ( |2.12[ ). This 
near saturation of the inequality is an indication of how well the constituent quark pic- 
ture works in this model. This is in agreement with the general argument concerning 
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the validity of the constituent quark picture in the large N c limit [17]. It is naturally of 
interest, then, to compare the constituent quark mass against the naive or the bare quark 
mass, which is done in fig. 4. In this figure, /i a a/2 is plotted against the bare quark mass. 
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fig. 4 The constituent quark masses of the fermion and the boson quarks computed from 
the meson masses compared against the naive quark mass, «/ /?/,&. "analytic" refers to the 
formula derived analytically from (\3.4) - 

The constituent quark mass is close neither to the bare quark mass m q nor the quantum 
corrected mass yjrriq — g 2 N c /7c, but is larger than m q by an amount of 0(1) when the 
mass is of the order of the QCD scale. This shows that even though the constituent quark 
picture seems valid, the constituent quark is not the naive quark but a "dressed" quark. 
When the quark mass is heavy, the difference between the constituent quark mass and the 
naive quark mass may be understood from analyzing the non-relativistic linear potential.0 



4. qq hadron vertex, confinement and asymptotic freedom 

In this section, we analyze some other physical properties of the "bf meson". The 
results are somewhat similar to the ff case || and the bb case || . In the bf case, we should 
note that unlike the ff and the bf cases, the "meson" is a fermion. 

First, let us analyze the qq scattering matrix. The equation satisfied by the qq scat- 
tering matrix is essentially that of the equation satisfied by the meson wave function ( |2.3|) . 

1 We would like to thank H. Kawai for pointing this out. 
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First we make use of the Lorentz structure of the scattering matrix and reduce it to its es- 
sential scalar component by defining T a ^{p,p' ;r) = (7_) a/ gT(p, p'; r). The reduced matrix 
element satisfies the equation 

T(p lP ';r) = ig 2 L(p_,p'_) + 2ig 2 N c J ' -^D(k)S(k - r)L(k-, P -)T(k,p'- : r) (4.1) 

Here we defined L(x, y) = (x + y)/(x — y) 2 and we denoted by S (p) the component of the 
fermion propagator that contributes to the scattering matrix r y-S(p)'j- = 2j-S(p). This 
equation may be solved in a manner closely related to that of the meson equation ( |2.3|) . 
Define </>(p_,p;r) as 

0(p_,p;r) = J dp + D(p)S(p-r)T(p 1 p';r) (4.2) 

so that 

/dk- 
— — L(k_,p_)(f)(k-,p;r) (4.3) 
Zn 

p + in Q4.2| ) may be integrated out explicitly to obtain the following equation for </> 



7T 2 

H(j)(x, y; r) = ii 2 ${x, y; r) + — L(x, y) (4.4) 

lM c r_x 

where H was defined in ( |2.3|) . We may obtain the solution to this equation using the meson 
wave functions as 

2 1 Z" 1 

(f)(x,x';r) = 2 / d y&k(x)Lp k (y')L(x',y') (4.5) 



Provided that {^fc} satisfies the following properties: 
1. The "Schrodinger equation" for the meson: 



H(p k {x) = ji 2 k (pk{x) (4.6) 



2. Completeness: 



3. Orthogonality: 



^2<Pk(x)<p k (y) = 8(x-y) (4.7) 



x 
k 



dxx(p k (x)(fi(x) = 5ki (4. 



o 



Corresponding properties in the ff case have been proven rigorously in [|18 
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Using the expression for the reduced scattering matrix (O) we may reconstruct the 
scattering matrix as 

T(x, x'; r) = % -^-L{x, x') + 2 $fc(aO$fc(s / ) (4.9) 

_ fc k 

where the g-g-meson vertex $>k( x ) is defined as 

a 2 {N \ 1/2 f 1 
® h{x) = l~ (if ) j Q dyL{x,y)(p k {y) (4.10) 

We see that the gg scattering matrix may be broken up into the g-g-meson vertex and the 
single particle poles corresponding to the mesons which are solutions to the bound-state 
equation. All the intermediate physical states to this order in 1/N C are mesons and we see 
that the colored particles are indeed confined. At higher orders in 1/N C , more complicated 
intermediate states such as two particle cuts will appear. 

To investigate the properties of the meson, we compute some matrix elements. Define 
a fermionic operator 

i 

where i is the color index. Denoting the fc-th meson state as hk, we may obtain some 
matrix elements such as 

/TV \ 1/2 r 1 

(0\F?\h k ) = -2i(-Z) dx(p k {x) 
N c \ rn f f <Pk(x) 



x 



Some correlation functions of these fermionic "currents" may also be obtained. Define 
the correlation function 

M^(g) = J d 2 xe^ x (o\TFf(x)FZ b (0)|o) (4.13) 

M__(g) = ^£ (jfWfcfr)) {AA ^ 

The behavior of this element in the deep inelastic region is 

2 

M__(g)~^ln^ for g 2 > m 2 , g 2 N c (4.15) 
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This agrees with the correlation function obtained in the free theory by computing the 
Feynman graph in fig. 5. This is, of course, none other than the statement of asymptotic 
freedom in this case. For comparison, the current matrix elements for the ff and the bb 

case are 



M—(q) 




(4.16) 



fig. 5 Correlation function of the fermionic currents in the free theory. 



5. Discussion 

In particle physics and in other areas of physics, such as condensed matter, we of- 
ten need to consider particles bound together by gauge interactions involving boson con- 
stituents. It is important to have a concrete example where such phenomenon is analyzed 
from first principles. In particular, the precise correspondence between QCD without 
quarks and string theory has been recently established in two dimensions []19| . It is of im- 



port to elucidate this correspondence when dynamical matter is coupled to QCD, putting 
the previous work on this subject [gj [2(J on a more rigorous footing. In this regard, it is 



crucial to understand in detail the properties of free strings, which are none other than 
the mesons in QCD. In this paper, we have analyzed the properties of bound states involv- 
ing boson constituents both analytically and numerically. We believe that this concrete 
model should contribute to the physical understanding of bound states involving boson 
constituents. 

There are a few intriguing aspects of our results which were not anticipated prior to 
computation: In mesons involving bosons, we found that it was not possible to construct 
massless bound states without additional interaction other than the gauge interaction. This 



is consistent with the result of |21[] where it is found that the phase transition between 
the broken symmetry phase and the symmetric phase in this model is of first order. This 
behavior differs from the ff case where zero mass quarks produced a massless meson. It 
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is possible that interactions may change some of these properties. Interactions that may 
be naturally added are the |</>| 4 interaction and the Yukawa interaction. There are, in 
addition, interaction terms that are renormalizable in 1+1 dimensions but not in 3+1 
dimensions, such as (if)i()) 2 and any gauge invariant scalar self interaction terms. 

In this work, we analyzed the QCD mass inequalities in 1+1 dimensions for mesons 
made from bosons and/or fermions. It is important to analyze this non-perturbative in- 
equality in a model when possible. We have shown how the inequality works quantitatively 
in 1+1 dimensions. An interesting aspect is that the QCD mass inequality is close to being 
saturated; as we have seen, to at most few percent relatively. Naively, there is no reason 
to expect that it should be so small and this is a sign that the constituent model works 
well even when the quark masses are light. We need to ask if these properties are artifacts 
of the model we used, namely the large iV c limit and the fact that the model is formulated 
in 1+1 dimensions. In the large iV c limit, quark loops are suppressed for group theoretical 
reasons and one may wonder if this is the cause for the saturation. First, even in two 
dimensional QCD in the large N c limit, the inequality is not exactly saturated and we 
know of no solid argument short of a concrete calculation that shows that the inequality 
is close to being saturated. Furthermore, in recent studies in 1+1 dimensions, it has been 
shown that even when the quark loops are not suppressed for group theoretical reasons, 
the lighter hadrons are very well approximated by a small definite number of partons, so 
that the quark loops are suppressed dynamically ||22|| . While it is unknown whether this 
feature is preserved in higher dimensions, we should keep in mind the successes of the 
constituent quark model JH]] |?3J . 

1+1 dimensional QCD has served particle physics well as a testing ground for various 
properties of QCD. In this paper, we have made an effort to present the method for 
computing the wave function and the spectrum of ff, bb and bf bound states in 1+1 
dimensional QCD in a concrete manner. We believe that this will be useful for further 
investigations in this field. 
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